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Abstract. — Here we present well-posedness results for first order stochastic differential inclusions, 
more precisely for sweeping process with a stochastic perturbation. These results are provided in 
combining both deterministic sweeping process theory (recently developed in [18] and [19]) and 
methods concerning the reflection of a Brownian motion ( [23) and [31[ ). In addition, we prove 
convergence results for a Euler scheme, discretizing theses stochastic differential inclusions. 

Resume. — Nous demontrons dans ce travail le caractere "bien-pose" d'inclusions differentielles 
stochastiques du premier ordre, plus precisement de processus de rafle avec une perturbation 
stochastique. Ces resultats sont issus de I'association de la theorie des processus de rafle 
deterministes (recemment developpee [18| et [19| ^ et de methodes concernant la reflexion d'un 
mouvement brownien ([23] et [31]). De plus, nous prouvons un resultat de convergence pour un 
schema d'Euler, discretisant ces inclusions differentielles stochastiques. 
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1. Introduction 

In this paper, we are interested in particular first order differential inclusions (namely sweeping 
process) with a stochastic perturbation. This work rests on the combining of two different 
theories: the first one about sweeping process and the second one about the reflection of a 
Brownian motion on a boundary. Let us first recall these two problems and related results. 
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Let B he a Banach space, X be a bounded time-interval, C : I ^ B he a set-valued map with 
nonempty closed values, and let / : X x ;B =^ be a perturbation. The associated sweeping 
process takes the form: 

' ^(i)+N(C(t),n(t))9/(t,u(t)) 

' u{t) G C{t) (1) 
, u{0) = uo , 

with an initial data uq £ C(0) and where N(C, x) denote the proximal normal cone of C at any 
point X. This differential inclusion can be thought as follows: the point u{t), submitted to the 
field f{t,u{t)), has to remain in the set C{t). 

This type of evolution problem has been extensively studied. It has been introduced by J.J. 
Moreau in 70's (see [29J) with convex sets C{t) of a Hilbert space and with no perturbation 
(/ = 0). To solve this problem, J.J. Moreau brings a new important idea in proposing a 
catching-up algorithm. 

Since then, some attempts have been made in the litterature to weaken the assumptions, for 
example to add a perturbation /, to weaken the convexity assumption of the sets, to obtain 
results in Banach spaces (and not only in Hilbert spaces). 

The perturbed problem in finite dimension {B = M'^) with convex sets C{t) (or complements of 
convex sets) was firstly studied by C. Castaing, T.X. Due Ha and M. Valadier in [11] . In this 
framework, they proved existence of solutions for ([1]) with a convex compact valued perturbation 
/ and a Lipschitzean set- valued map C. Then in [12j, C. Castaing and M.D.P. Monteiro Marques 
considered similar problems in assuming the upper semicontinuity of / and a "linear compact 
growth" : 

f{t,x) C {3{t){l + \x\)B{0,l) , y{t,x) € / X M'^. (2) 

Moreover the set-valued map C was supposed to be Hausdorff continuous and satisfying an 
"interior ball condition": 

3r > 0, 5(0, r) C C{t) , G /. (3) 



Later the main concept which weakens the convexity property of sets C{t), is the notion of 
"uniform prox-regularity" . A set C is said to be uniformly prox-regular with constant tj or r]- 
prox-regular if the projection onto C is single- valued and continuous at any point distant at most 
r] from C. 

The uniform prox-regularity assumption was made in numerous works dealing with sweeping 
process. The Hilbertian case without perturbation (/ = 0) was firstly treated by G. Colombo, 
V.V. Goncharov in [16| . by H. Benabdellah in [2] and later by L. Thibault in |37) and by G. 
Colombo, M.D.P. Monteiro Marques in [17]. In j37j. the considered problem is 

-du G N(C(t),u(t)) 

, ' (4) 

u{To) = Uo , 

where du is the differential measure of u. The well-posedness of dH is proved under the same 
assumptions as previously excepted ([3]). 

In an infinite dimensional Hilbert space B {B = H), the perturbed problem is studied by 
M. Bounkhel, J.F. Edmond and L. Thibault in fH [371 [18l [19] (see Theorem [SJ]) and recently 
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by the authors in a Banach space in [7J. For example in [19], the well-posedness of 

-du G N(C(t),u(t)) + f(t,u(t))dt 

(5) 

m(0) = uo , 

is proved with a set-valued map C taking rj-piox regular values (for some rj > 0) such that 

\dc{t){y) - dc(s){y)\ < K]s,t]) , yyeH,ys,tei,s<t (6) 

where is a nonnegative measure satisfying 

sup//({s}) < ^. (7) 

S&I ^ 

All the proofs rest on the algorithm developed by J.J. Moreau with additional arguments to deal 
with the prox-regularity assumption. 

We now want to recall results about the reflected Brownian motion in a set. Before that, we 
refer the reader to jS], IIP) , where C. Castaing proved some results about existence of solutions 
for sweeping process with a convex set C{t) which is stochastically perturbated. Here we want 
to add a stochastic perturbation in the differential inclusion without changing the deterministic 
time-evolution of C. 

Reflected Brownian motion 

We consider a closed set C in and we look for solving stochastic differential inclusions, 
describing the time-evolution of a Brownian motion (inside C) with a reflecting boundary dC . 
Let {Bt)t>o be a M"^- valued Brownian motion, then the path is given by a stochastic process X 
involving the following stochastic differential inclusion: 

dXt+^^{C,Xt) 3 dBt 

XteC 

Xo = Uq, 

where uq is the starting point. One of the difficulty is to give a precise sense to "N(C, X^)". 

The first well-posedness results have been investigated by many authors (see for example 
A.V. Skorohod [33l [34], N. Ikeda and S. Watanabe [21l H] and N. El Karoui [22] ...) for 
C a half-plane (or a half-line). Then the problem was treated for a smooth set C by D.W. 
Stroock and S.R.S. Varadhan [35] and later by A. Bensoussan and J.L. Lions [3]. H. Tanaka 
succeeded in getting around the smoothness assumption in the case of a convex set in [36] . In 
[23] ■ P.L. Lions and A.S. Sznitman gave the first result for bounded uniformly prox-regular sets 
(without this terminology) in assuming an extra assumption of "admissibility" . Few years later, 
Y. Saisho has extended the proof for unbounded prox-regular sets in [31j . They considered 
the associated deterministic Skorohod problem, which consists for a continuous function h on 
a time-interval I in finding a pair of continuous functions (x, k) defined on I with k G BV{X) 
satisfying: 

Vt G X, x{t) + k{t) = h{t) (8) 
and ^ ^ 

\m)= ! Us)^dcd\k\{s), k{t)= [ Csd\k\is), (9) 
Jo Jo 

with G N(C,x(s)). 

The equation ([8]) corresponds to an integral version of ([1]) with a constant set C{t) := C and 
h{t) = Jq f{s,u{s))ds. Indeed ([9]) specifies that the support of the differential measure dk 
is the set of the snapshots t when x reaches the boundary dC and gives a precise sense to 
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"d/c G N(C, . So the Skorohod problem can be thought as an integral version of a sweeping 
process. Using this deterministic problem, the previously cited works deal with the following 
stochastic differential inclusion 

dx{t) + N{C,x{t)) 3 f{t,x{t))dt + a{t,x{t))dBt 

< u{t) G C 

. u(0) = no , 

where / and a are some perturbations and {Bt)t>o is a real Brownian motion. This problem 
can also be seen as a stochastic perturbation of the "sweeping process" by a constant set C. 
We would like to finish this brief state of art by referring the reader to |32] for a work about 
reflecting Brownian motion in a set corresponding to the complement of a collection of balls. 
We are specially interested by this example and we are looking for extending this result to more 
general situations with time-dependent constraints. 

Associated numerical schemes 

Concerning the deterministic sweeping process, the existence results are obtained by the 
convergence of the so-called "catching- up algorithm" . More precisely, in considering some sub- 
division {Jk)k of the time-interval, the set-valued map C is approached by a piecewise constant 
multifunction taking value on = [tk-,tk+i[. The associated discretized solution u is defined 
as follows: 

Vt G Jfc+i , u{t) = Uk+i = Pcfe+i [uk + {tk+i - tk)f{tk,Uk)] , 

with no fixed to the initial value. Then it is proved that under the above assumptions, the 
sequence of discretized solutions is of Cauchy type and also converges to a function, which is 
the solution of the continuous problem ([T|). 

About the first order stochastic differential inclusions, we would like to present the work of E. 
Cepa [13] and F. Bernardin They consider equations taking the form 

' dx{t) + A{x(t))dt 3 f{t, x{t))dt + a{t, x{t))dBt 

< u{t) G C 

^ n(0) = no , 

where A is a time-independent maximal monotone operator. They define Euler numerical scheme 
(using the resolvents and Yosida's approximation of a monotone operator) and prove its conver- 
gence. In our case, (even for time-independent set) we are interested in uniformly prox-regular 
set C and so, it is well-known that the associated proximal normal cone (as multivalued operator) 
could be not maximal monotone. It only satisfies a weaker property of hypomonotonicity. 

Framework for Stochastic perturbation of sweeping process and main results 

We now come to our contribution. As explained before, the theories of sweeping process and 
reflected Brownian motion are based on the uniform prox-regularity property of the sets. In 
order to combine these two theories, we consider a set-valued map C(-) which is admissible and 
regular (in extending these defintions to set- valued maps, see Definitions 13.11 and 13. 3|) . 

We begin this work by describing some abstract results about uniformly prox-regular sets in 
order to understand our required assumption. Then, we propose general results about a set- 
valued map C. Extending the assumptions of |23) to the time-dependent framework, we prove 
the two following results in Section [3l 
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Theorem 1.1. — Let X he a hounded time interval and C : X ^ M admissihle, regular set- 
valued map varying in an ahsolutely continuous way (see Definitions \3.1l \3.3\ and \3.4\ )- 
Let /, fj : I X M"^ — )• M'^ be bounded and Lipschitz with respect to the second variable. Then the 
following problem 

' dx{t) + N(C(t), 3 f{t, x{t))dt + cr(t, x{t))dBt 

< u{t) G C{t) (10) 
, u(0) = uo, 

(given by a M-Brownian motion {Bt)teX on a probability space P)j is well-posed in 

L^(r2, L°°(T)). That means: for every initial data uq G C(0), there exists one and only one 
process {Xt)t^x £ L'^{Q.,L°°(1)) solution of \21\) (in the sense of pathwise uniqueness). 

Moreover we will obtain stability results, which permit to understand the behavior of the stochas- 
tic solutions when the perturbation a tends to (see Theorem 13. lip . 

If the regularity of C(-) is not assumed, well-posedness results still hold (but not necessary in 
the space L^{n, L°°{I))): 

Theorem 1.2. — Let X be a bounded time interval and C :X ^ M*^ admissible set-valued map 
varying in an absolutely continuous way (see Definitions \3.1\ and \3.4\ )- 

Under the same assumptions on f and a, the problem Iil0\) has a unique solution (in the sense 
of pathwise uniqueness). 

In Section U we prove convergence results for a Euler scheme, discretizing the stochastic differ- 
ential inclusion ()10p . 

In Section [5l we check that the main assumptions are not too strong in applying these results 
to particular set- valued maps C(-) given as the intersection of complements of convex sets. 
Moreover an application to a model of crowd motion is described. 

This work was briefly presented in a proceeding [6]. 



2. Preliminaries about Prox-regular sets 

We emphasize that the different notions defined in this section can be extended in the case 
of an infinite dimensional Hilbert space H. Here we only deal with a Euclidean framework. We 
denote by B the unit closed ball. For a subset C of M*^, we write dc for the distance function to 
this set: 

dc(x) := inf \y — x\. 

Definition 2.1. — Let C be a closed subset ofW^. The set-valued projection operator Pc is 
defined on by 

Vx G Pc{x) := {y G C, \x - y\ = dc{x)} . 

Definition 2.2. — Let C be a closed subset ofW^ and x G C, we write N(C, x) for the proximal 
normal cone of C at x, defined by: 

N(C,x) := {v G M'^, 3s > 0, x G Pc{x + sv)\ . 
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We now come to the main notion of uniformly prox-regular set. It was initially introduced by 
H. Federer (in |20j ) in finite dimensional spaces under the name of "positively reached set". 
Then it was extended in infinite dimensional space and studied by F.H. Clarke, R.J. Stern and 
P.R. Wolenski in [15] and by R.A. Pohquin, R.T. Rockafellar and L. Thibault in }30l . 

Definition 2.3. — Let C be a closed subset ofM.'^ and i] > 0. The set C is said r] -prox-regular 
if for allxeC and v € N(C, x) \ {0} 

B + r/|^,r/^ n C = 0. 

Equivalently, C is r] -prox-regular if for all y ^ C , x ^ dC and v € N(C, x) 

(y -x,v) < |^|x - yp. (11) 

Remark 2.4- — We refer the reader to jl5|, [14] for other equivalent definitions related to the 
limiting normal cone. Moreover we can define the notion using the smoothness of the function 
distance dc (see [3jQjJ. 

This definition is very geometric, it describes the fact that we can continuously roll an external 
ball of radius smaller than rj on the whole boundary of the set C. The main property is the 
following one: for a r/-prox-regular set C and every x satisfying dc{x) < rj, the projection of x 
onto C is well-defined (i.e. Pcix) is a singleton) and the projection is continuous. 

Definition 2.5. — Let A : ^ be a set-valued operator on M.'^. We write D{A) for its 
domain: 

D{A) := |x G M.'^, A{x) / 0} 

and T{A) for its graph: 

r{A) := {(x,a) G M^"^, x G D{A), a G ^(x)} . 

The operator A is said "hypomonotone" if there exists a constant 5 > such that for all x,y (z 
D{A) and all (q,/3) G A{x) x A{y), we have 

{a-(3,x-y) > -S[\a\ + m\x-y\\ (12) 

Such an operator is called hypomonotone with constant 6 and it is maximal if for all x,q G M'^ 
then 

(Vy G D{A), V/3 G A{y), {a - f],x - y) > -S [\a\ + |/3|] |x - y\^) ^ a G A{x). 

Then we refer the reader to the work |30| of R.A. Poliquin, R.T. Rockafellar and L. Thibault 
for the following result: 

Proposition 2.1. — Let C be a closed subset. Then C is r] -prox-regular if and only if the 
proximal normal cone N(C, •) is a maximal hypomonotone operator with constant Moreover 
L»(N(C,-)) = C. 

The following definition comes from the work of A. A. Vladimirov (Section 3 in [40) ). We extend 
it to hypomonotone set-valued operators. 

Definition 2.6. — For A and B two hypomonotone operators with constant 5, we define 

, / . {a- (3,y-x) -25{\a\ + |/3|) |x - 

dv{A,B):= sup sup , , , , , • 

(x,a)Gr(A) (s/,/3)Gr{B) i + |a| + \p\ 
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Note that dy is not a distance as the triangle inequahty is not satisfied. However, we have the 
following inequality: 



Proposition 2.2. — Let Ci and C2 be two rj-prox-regular sets o/M and write du the Hausdorff 
distance 

dH{Ci,C2) ■■= sup |dci(y) -dc2{y) \ ■ 



Then with S = we have 

dv(N(C7i, •),N(C2, •)) < dniCi, C2) + -dH{Ci,C2f. 

Proof: The proof is the same one as for Lemma 3.4 in [40] (dealing with convex sets). We 
detail it for an easy reference. Let (x,a) G r(N(Ci, •)) and G r(N(C2, •)). We know that 

there exists a point x G Pc^ (y) G Ci such that |y — i| < d/f(Ci, C2). Then the hypomonotonicity 
property of the proximal normal cone implies that 

{a,x — x) < — lallx — 

27] 

1 2 

- 2?]'"' ^'^ ~ 2/1 + |y - x\] 

< -\a\ [\x - + ly -xp] , 

hence 

(a, y - x) < |a|(i_f/(Ci, C2) + -|a||x - + -laldj^fCi, C2)^. 

Since an analoguous inequality is satisfied by /3, the desired result is obtained by summing both 
inequalities. □ 

Corollary 2. 7. — Let C : X ^ he a set-valued map taking r]-prox-regular values. If t ^ 
C{t) is absolutely continuous (for the Hausdorff distance), then the operator-valued map t — > 
N(C(t),-) is absolutely continuous for dy ■ 

According to the work of A. A. Vladimirov ( |40j ). an absolute continuity of the operator-valued 
map t N(C(t), •) is the appropriate assumption to solve the differential inclusion: 

^{t)+NiC{t\x{t))3fit,xit)). 

That is why, the previous proposition suggests to suppose an absolute continuity for the map 
t ^ C{t) (which is exactly the same assumption as done in [18t, il9j ) . 

We finish this section with the following result: 

Proposition 2.3. — Let C be a rj-prox-regular set in W^. Then for all e < r]/8, the set 

C, := |x G M"', dc{x) < e} = C7 + eB 

is uniformly r]/8-prox-regular. Moreover the distance function: dc^ is in a neighborhood of 
dC^. So the boundary dC^ is a -manifold. 

Remark 2.8. — This proposition is based on general results due to F. Bernard, L. Thibault 
and S. Zlateva in [4] and so it can be extended in a uniformly convex Banach space framework. 
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Proof: We will use the notations of [4j, mainly for a closed set C, we define 

Ec{l) := G R"^, dc(a;) > ^} • 

^Prom Theorem 6.2 and Lemma 6.9 of [5, we deduce that for I = rj/4 the set Ec{l) is Z-prox- 
regular. Then Lemma 12.91 (later proved) ensures that for all e G (0, //2) we have 

C, = C + eM = EEaii){l-e). 

Analogously, we deduce that the set is (/ — e) prox-regular and so is ?7/8-prox-regular. 
Moreover it can be checked that these sets have a boundary. Indeed Lemma 12.91 vields 



G M°', dc{x) - e = o| . 



Applying Theorem 6.2 of we deduce that dc is in a neighbourhood of dC^ and so dC^ is 
a C"^ manifold. □ 

The proof is ended provided that we prove the following lemma: 

Lemma 2.9. — Let C he a -q -prox-regular set and I = r//4. Then for all e G (0, 1/2) 

:= [x G dc{x) < e} = Eec(i){1 - e). (13) 

Moreover 

dC, := {x G M'^, dc{x) = e} . (14) 



Proof: We first check the two embeddings of (|13p . First as the function dc is 1-Lipschitz then 
for e < Z, it yields 

Vx G C + eB, Vu G ^c(0> ^ - e < dc(M) - dc{x) <\u- x\. 

So we deduce that dE(^{i){x) > / — e for all a; G C + eB, which gives 

Ce = C + eBci?^^(i)(/-e). 

For the reverse inclusion, we have to use the prox-regularity assumption. Let us take x G 
Ee(j{i){1 — e). We may suppose that x does not belong to C, else x G C C C + eB. So let us 
consider xq G Pc(x). Then a; — xq is a proximal normal vector at xq and so we know that 

X — Xn ^ 

u ■.= Xo + l-. ^ G Ecil). 

|x — Xol 

Here we have used that / < rj/2 to get dc{u) = I. Then as x G EE(y(i){l — e), we deduce that 

^ — e < |n — x| = |/ — |x — Xol I = / — |x — xo|, 

where the last equality comes from the fact that x does not belong to Ec{l)- The inequality 
yields |x — xo| < e, which ends the proof of 

C + eMZiEEa{i){l-e). 

Finally, it remains us to prove ()14p . First thanks to the continuity of the distance function, it is 
obvious that 



dC, C {x G M'^, dc{x) = e} . 



The previous reasoning (based on the prox-regularity of the set C) yields the reverse embedding 
and also concludes the proof. □ 
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3. Well-posedness and Stability results for stochastic sweeping process 

In all this section, we consider a set-valued map C. We first define the required assumptions, 
under which well-posedness results for stochastic sweeping process can be proved. We want to 
study the stochastic differential inclusion (jlOp which can be seen as the equation of a reflected 
Brownian motion onto the moving set C(-). In order to apply the results of P. L. Lions and A. 
S. Sznitman [ 23] , the sets C{t) are supposed to satisfy some properties. 

We refer the reader to [23] for the following definitions without considering the time variable. 
Here we add the time-dependence. 

Definition 3.1. — The set-valued map C is said admissible on [0,T] if it takes uniformly 
prox-regular values (with a same constant) and if there exist 5,r,T > 0, and for all t G [0, T] 
sequences {xp)p and {up)p with \up\ = 1 and Xp G C{t) such that for alls € [0,T] with \t — s\ < t, 
{B(xp,r))p is a bounded covering of the boundary dC{s) and 

yp, yy £ dC{s)nB{xp,2r), yv eN{C{s),y), {v,Up)>d\v\. (15) 

Remark 3.2. — The original definition in |23] makes appear an extra assumption: there exists 
a sequence of approaching "smooth sets" satisfying a uniform bound of prox-regularity. As 
explained in |31j . the existence of such approaching smooth sets is not really necessary to prove 
the solvability of the Skorohod Problem. Moreover, due to the recent works about prox-regular 
sets, we know that such an approaching sequence always exists in a very general framework (see 
Proposition \2. 3\) . 

The second important property (used in [23]) is the following one: 

Definition 3.3. — A set-valued map C is said regular on I if there exists a function $ G 
Cl{I X W^) satisfying for all t G [0,r], x G dC{t) and v G ^{C{t),x) 

{V^^{t,x),v) < -ii\v\, 

with > 0. 

According to Remark 3.1 of j23j and Lemma 5.3 of [3,1], we know that an admissible set is 
locally regular. Moreover, Y. Saisho has proved in |31] that the local regularity is sufficient to 
obtain well-posedness results. 

Definition 3.4 ([181 119] ). — A set-valued map C{-) is said to vary in an absolutely continuous 
way, if there exists an absolutely continuous function v such that for all t,s (zl 

dH{C{t),C{s)) < \vit)-vis)\, 

where dn is the Hausdorff distance. 

3.1. Deterministic Skorohod problem: an extension of sweeping process. — Before 
solving the stochastic differential inclusion ([10]), we study the associated Skorohod problem. Let 
h :I he a continuous function satisfying /i(0) G C(0). We say that a couple of continuous 

functions (x, k) on Z is a solution of the Skorohod problem (SkP,/i) if: 

• for all t G I, x{t) G C{t) 

• the function k is continuous and have a bounded variation on I 

• the differential measure dk is supported on {t G I, x{t) G dC{t)}: 

\k\{t) = I l^(^s)&dC(s)dms), k{t) = f as)d\k\{s), (16) 
Jo Jo 

with C(s) G N(C7(s),x(s)) 
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• for all t € X, we have 

x{t) + k{t) = h{t). (17) 

Here we denote by \k\(t) the total variation of the function A: on [0,t]. By extension, a continuous 
function x is said to be a solution of (SkP,/i) if there exists a function k such that (x, k) satisfies 
the previous properties. 

This subsection is devoted to the study of the deterministic problem (SkP,/i) (defined by (fT6|) 
and ()17p ). By following the ideas of |23^ I31j . we also begin with the following proposition: 

Proposition 3.1. — Consider an admissible set-valued map C varying in an absolutely con- 
tinuous way and assume that for all h G C°°{I,M.'^) with h{0) € C(0), there exists a solution 
{x,k) to the Skorohod problem (SkP,h). 

Then for all h G C^[X,W^) with h{0) € C{0), there exists a unique solution {x,k) to the Sko- 
rohod problem (SkP,h). Furthermore the mapping [h — > x) from C^{X,W^) into itself is Holder 
continuous of order ^ on compact sets. 

Proof: We refer the reader to Theorem 1.1 of |23j for a detailed proof of such proposition in the 
case of a constant set C. It is based on the hypomonotonicity and the admissibility properties. 
First step: Uniqueness. 

First the uniqueness is "as usual" a direct consequence of the hypomonotonicity property and 
Gronwall's Lemma. We have to be careful as we are working with a function k which is only 
assumed to have a bounded variation. 

We recall the following version of Gronwall's Lemma (due to R. Bellman [T]): 

Lemma 3.5. — LetX := [0, T] be a closed bounded interval. Let u be a non-negative measurable 
function defined on Z and let fi be a locally finite non-negative measure on I. Assume that 
u € L^(Z, dfi) and that for all t ^ I 

u{t) < / u{s)dfi{s). 

J[o,t] 

If the function t — )• /x([0, t]) is continuous on Z, then for allt^Z 

u{t) = 0. 

Let us take two solutions x and x (we associate the corresponding functions k and k). We study 
the error term: 

z{t) := x{t) - x(t) e{t) := |z(t)p. 

The difference function z is also solution of the following differential equation (in the sense of 
time- measure) : 

dz{t) = -dk{t) + dk{t). 

As k and k have a bounded variation, it comes 

de{t) = 2{z{t),dz{t)) = 2{z{t), -dk{t) + dk{t)). 

^From the hypomonotonicity of the proximal normal cones N(C(t),-) (due to the r/-prox- 
regularity property of the sets C (t) , see Proposition 12. ip , we deduce that 

deit)<^{\dk\it) + \dk\it)y 

Then as k and k are assumed to have finite variation on Z and e(0) = 0, we deduce that e = 
thanks to Lemma [331 (with fj. = \dk\ + \dk\) . The proof of uniqueness is also concluded. 
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Second step: Existence. 

Let h G C°°(I,R'^) and {x,k) a solution to the Skorohod Problem (SkP,/i). We follow the 
reasoning, used in the proof of Lemmas 1.1 and 1.2 of [23j, in pointing out the difficulties raised 
by the time-dependence of C. 

First, Lemma 1.1 (6) of [23] still holds as it only rests on the uniform prox-regularity of sets 
C{t). We detail how Lemma 1.1 (7) of [23j should be modified to take into account the time- 
dependence of C. Let r, r be given by the admissibility property (Definition 13. ip . On [0, r], we 
denote by {Oi)i the sets Oi := B{xi, 2r) n C(0) and Oq the following set 



Oo = CiO) \ 



|j5(x„r) 



We let Ti := inf{t G [0,r], x{t) ^ OiJ, where io is such that h{0) = x(0) G Oi,,. Then either 
x(Ti) G Oq and we set ii = 0, or x{Ti) G B{xi^,r) for some ii. In this way we construct by 
induction, im and Tj^ such that if Tm < r, x{Tm) G B{xi^,r) or x{Tm) G Oq, T^+i := inf{t G 
]Tm,r], x{t) ^ Oi^}. Using the admissibility properties, we obtain such estimations for the 
variation of k: let s,t G [rm-,Tm+i]-, s < t for some m) 

\k\{t)-\k\{s)<-^{\h{t)-h{s)\ + \x{t)-x{s)\). 

Moreover if = then for all t G \TrmTm+i\ x{t) ^ dC{t) and the variation is constant 
on [TmjT'm+i] else for all t G [Tm,Tm+i] x{t) G B{xi^^,2r) and \x{t) — x{s)\ < 2r. We finally 
conclude that 

\k\{Tm+l)-\k\{Tm)<K, (18) 

where X is a numerical constant depending on \h{Tm+i) — h{Tm)\- Then we can repeat the 
proof of Lemma 1.2 in |23j in order to give an upper bound to the variation |A;|(r). However 
the time-dependence of C makes appear a new quantity in (7) of |23j . More precisely, we get 
for s,tG [T^,T^+i]^ 

h{s)\ 



\<t) - x(s)|2 < sup \h{-) - h{s)\^ + sup \h{- 

[s,t] [s,t] 



+ i\k\it) - \k\is)) 



supdHiC{-),Cis))+snpdH{Ci-),C{s)f 

[s,t] [s,t] 



(19) 



Indeed, in the proof, we have to estimate for u G [s, t] 

2{x{u) — x{s),dk{u)) — rid\k\{u)\x{u) — x(s)p. 

To use the hypomonotonicity of the proximal normal cone (see (jlip l. we have to write x{s) = 
y{u) + z{u) with ?/(?x) G Pc{u){x{s)). As x(s) G C(s), \z{u)\ < dH{C{u),C{s)) and this operation 
makes appear a rest which is bounded by dH{C{u),C{s))d\k\{u). 

The two last terms in ()19p can be assumed as small as we want since C varies in an absolutely 
continuous way and the variation of k is bounded by (jlSp . If r is small enough (with respect to 
r), then by definition of T^+i and we deduce that for some constant c 

< |x(T^+i) -x(T^)|2 < c sup \h{-) - h{s)f + c sup \h{-) - h{s)\ + r'^ /2 

lTm,Tm+l] [Tm,Tm + l] 

which gives 

r^< sup \h{-) - h{s)\^ + sup \h{-)-h{s)\. 

lTm,Tm + l] lTm,Tm+l] 



For two quantities A, B, we write A < B if there exists a constant C such that A < CB. 
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As done in [23j . we can also conclude that the length \T„i+i — 7m | is not too small, according 
to the uniform continuity modulus of h and so that the collection of indices m is finite and 
its cardinal is bounded by a constant depending on h. Finally, the total variation |A:|[o,r] •= 
\k\{T) — \k\{0) = |A;|(r) is controlled by \\h\\oo and the uniform continuity modulus of h. 
The same reasoning holds on [r, 2r]. Let us remark that if {x,k) is the solution of (SkP,/i) on 
[0, 2r] then (x, k — k{T)) is the solution of (SkP,/i — A;(r)) on [r, 2r]. So we know that the variation 
\k — /c(r)|[^2T] = l^l[r,2r] Can be bounded by the norm of /i — A;(r) (which is controlled by 
the L°° norm and the continuity modulus of h) and its uniform continuity modulus (which 
is the same as for h). Consequently, the total variation |fc|[o^2r] is bounded with the help of 
h G C'^([0, 2t]). By iterating (at most T/r steps), we conclude that for h G C°°, the (assumed) 
solution (x, k) satisfies 

lkllL°°([0,T]) + l^l[0,T] < 

where K only depends on the L°° norm and the uniform continuity modulus of h. This key-point 
permits us to take the limit and to prove Proposition 13.11 by density arguments as in |23j (the 
proof can be easily adapted with a time-dependent set C). □ 

Theorem 3.6. — Consider an admissible set-valued map C(-), varying in an absolutely con- 
tinuous way. For all h G C^(X,W^) with h{0) G C(0), there exists a unique solution {x,k) to 
the Skorohod problem (SkP,h). Furthermore the mapping {h — )• x) from C^{Z,W^) into itself is 
Holder continuous of order ^ on compact sets. In addition if h ^ BV{I) then x G BViX) and 

d\k\t<d\h\t + d\v\u (20) 

where v is given by Definition \3.4\ 



Proof: According to the previous proposition, we just have to deal with smooth functions h. 
Indeed for h G VF^'^(X), the Skorohod problem (SkP,/i) has already been treated since it corre- 
sponds to the so-called sweeping process. We refer the reader to |18|, 119] for more details about 
these differential inclusions. We deal with Lipschitz moving sets C{t). So applying Theorem 1 of 
|18j . we know that there exists a pair of function (x, k) G VF^'^(X)^ such that almost everywhere 
on I 

dx dk dh 
dt dt dt 

and for almost every t G I 

dh 

-{t)eN{Cit),x{t)). 

Moreover, ([20]) holds. 

This gives us a solution to the Skorohod Problem (SkP,/i) for every smooth function h G W^'^{Z). 
Then we conclude the proof of the theorem, thanks to Proposition 13.11 □ 

3.2. Sweeping process with a stochastic perturbation. — In this section, we consider 
the Euclidean space M.'^ (equipped with its Euclidean structure), a probability space {ft,T,F) 
endowed with a standard filtration {J-t)t>o and a standard M-valued Brownian motion {Bt)t>o 
associated to this filtration. We denote by E(X) the expectation of a random variable X, 
according to this probability space. 

We fix a bounded time-interval X = [0, T] and denote by B the Banach space of time-continuous 
J't-adapted process X satisfying 

r 

\\X\\b:=^ suplXfl"^ <oo. 
Itex 
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Let / : I X M'^ — > M°' and a : X x M"' ^ M'^ be two maps L^. We look for solving the following 
stochastic differential inclusion on X: 

( dXt + ^{C{t),Xt) 3 fit, Xt)dt + a{t, Xt)dBt 

(21) 

[ Xq = Uo, 

where uq € C(0) is a non-stochastic initial data. Let us first give a more precise sense to this 
differential stochastic inclusion. 

Definition 3.7. — A continuous process {Xt)tei is a solution of \21\) if there exists another 
process [Kt)t<^x such that: 

a) {Xt)tel "is a -valued process taking values in C{t) and Tt-o-dapted with continuous sample 
paths; 

b) {Kt)t^x is a W^-valued process, Tt-adapted whose sample paths are continuous and have a 
hounded variation on X; 

c) the following stochastic differential equation is satisfied 

dXt + dKt = fit, Xt)dt + ait, Xt)dBt; (22) 

d) the initial condition is verified: Xq = uq P-a.e.; 

e) the process dKt is supported on {t, Xt G (?C(t)}.- 

\K\t= f lx.^dC{s)d\K\s, Kt= f iis)d\K\s, (23) 
Jo Jo 

with Cis) G N(C7(s),X,). 
The point e) gives a precise sense to "dK^ G N(C(t),X()". 

Using the stochastic integral, we can rewrite (j22|) as follows: we are looking for processes iXt)t(^x 
and iKt)tex satisfying that there exists a measurable set Qq C O of full measure such that for 
all u; G r^o and all t G X we have: 



Xt = uo+ f fis, Xs)ds + / ais, Xs)dBs - Kf 
Jo Jo 



(24) 



3.2.1. Well-posedness results for \21\) . — Now we come to our main results. We first assume that 
the set-valued map C is admissible and regular in order to follow the ideas of P.L. Lions and 
A.S. Snitman in [23j and obtain well-posedness results in the space L^(r2, L°°(Z)). We then 
describe results without requiring the regularity assumption following the ideas of Y. Saisho in 



Theorem 3.8. — Consider an admissible and regular set-valued map C, varying in an abso- 
lutely continuous way. 

Let f ,a : X xW^ ^ be bounded and Lipschitz with respect to the second variable: there exists 
a constant L such that for all t ^ X and x, y G M'^ 

\ait,x)-ait,y)\ + \fit,x) - fit,y)\ < L\x - y\ 

and 

\fit,x)\ + \ait,x)\ <L. 



'^'We have chosen an R'^-valued function a with a real Brownian motion (-Bt)t>o. Indeed all the results and the 
proofs hold for a real function a with an R'*-valued Brownian motion. 
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Then \21\) is well-posed in L^(i7, L°^{Z)). That means: for all initial data uq G C(0), there exists 
one and only one process {Xt)tei € L^(r2, L°°(X)) solution of f2l\) (in the sense of pathwise 
uniqueness) . 

Proof of Theorem 13. 8t For all X S ^B, we write F{X) the unique solution of the Skorohod 
Problem (SkP, h) with 

F{X)t + Kt = h{t) ■.= UQ+ f\{s,Xs)dBs+ I f{s,X,)ds. 

Jo Jo 

Due to Theorem 13.61 and the continuity of the stochastic integral, the unique solution F{X) 
exists. This well-defined map F satisfies some properties as pointed out by the next proposition. 

Proposition 3.2. — The map F takes values in B. There exists a constant k such that for all 
processes X, X' £ B, 

rT 



F{X)-F{X')l<k 



E 



\X -X 



/||4 



lL-([0,t]) 



dt. 



(25) 



The constant k only depends on \Z\ = T and on the above constants about C(-), / and a. 

Consequently, for all integer p > 2, we have (by writing F^^^ = F o .. o F for the iterated 
function) : 



F^P\X) - F^'P\X') 



< k 



< k' 



E 



||i7(p-i)(x) -F(p-^)(X') 



Il°°([oa]) 



dti 



T ft 



E 



Jo 
T rti 



Wf'^p-^Xx) - f'-p-^^x') 



lL°=([0,t2]) 



dt2dti 



< kP 

TPkP 

< 



Jo 



tp-l 



E 



\X - ^'||ioo([o,tj,]) 



dtp ■ ■ ■ dti 



p 



\X -X 



/||4 

Is- 



We also deduce that for a large enough integer p (only dependent on the above constants), the 
map F^P'^ is a contraction on B. By the well-known fixed-point Theorem, we also know that F 
admits one and only one fixed point X £ B. It is obvious that X is then a solution of ()24p and 
()2ip . So we conclude the proof of Theorem 13.81 provided that we prove Proposition 13.21 □ 

Proof of Proposition 13. 2t 

We only deal with ()25p . a very similar reasoning permits to check that F takes values in B. 
We follow the ideas and the estimates of |23j . Let X,X' G B, we denote Y := F{X), Y' := 
F{X'), K,K' the associated processes and $ the function given by Definition 13.31 Then Ito's 
formula implies 



^t, Yt) = <^{0, no) + ^{s, Ys)ds + J\^x^s, Y^), /(s, Xs))ds 



{a{s,Xs),S/,<^>{s,Ys))dBs 



{S/,^{s,Y,),dK,) 



+ 



{a{s,Y,),H,cl){s,Ys)a{s,Ys))ds. 



(26) 
(27) 



We denote the Hessian matrix by Hx. So with a := — it comes 

ri 

expi=^mt,Yt) + <^{tX))\ \Yt-Y;\^ = 1 + 11 + 111 
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with 



I :=2 exp — (c^>(s,y,) + <I>(s,yi)) {Y, -Y^,a{s,X,) - ais,Xi))dB. 



a 



(28) 



+ {Ys - y;, f{s, Xs)ds - f{s, X',)ds - dKs + dKi) 

+ 1^ exp Ys) + y/)) I tr [{a{s, X^) - a{s, X'M^{s, X,) - a{s, X',)Y] ds 

corresponding to the differentiation of the square quantity \Yt — If'P, 



a 



exp —ms, Ys) + Hs, y/)) |y, - y,r (v,$(s, y,), f{s, Xs))ds 



+ {V^<!>{sX)J{s,X',))ds + {a{s,Xs),V.^s,Ys))dBs + {a{s, X',),V^^sX))dBs 
^{a{s, Xs),H^^{s, Y,)a{s, Xs))ds + ^{a{s, K),H,4>{s, Y^)a{s, X',))ds 
- {V^<!>{s,Ys),dKs) - {V^<!>{sX),dK) 

+ —(s,Y,)ds + —{sX)ds 



1 

2^ 



-1 



exp —iHs, Ys) + Hs, y;)) \Ys - y. 



a 



rl\2 



(29) 



(v,«>(s, y,)a(s, Xs) + v,$(s, y;)a(s, X0)2 



ds 



corresponding to the differentiation of the exponential quantity (with the Ito's additional term 
in dllD) and 



in 



a 







-1 



exp — (ci>(s,y,) + cD(5,y/)) 



a 



(30) 



(y, - y;, a(s, X,) - <j{s, X',)) [VMs, Ys)a{s, X^) + V,.$(s, y;)a(s, X^)] ds 



corresponding to the variational quadratic term due to the product between the square quantity 
and the exponential one. Then we recall that $ G C|, by Definition 13.31 and that C{t) is 
ry-prox-regular, we have 

1 



a 
1 



-(v,$(s, y,), dK,)\Ys - y,V - {Ys - y;, di^.) < 



-{vMsX).dK)\Ys - y,f - (y. - y;,di^:) < o 

a 

in the sense of a nonnegative time-measure. So from the boundedness and the Lipschitz regu- 
larity of a and /, we deduce that 



\Yt-Yl\^< 



-1 



a 



exp — ($(s, n) + $(s, y;)) (y, - y;, ct(s, X,) - a(s, K))dB, 



f Ws\Ys - [(a(s, X,), V,.cl>(s, y,)) + (a(s, X^, V,.$(s, y/))] di?, 

JO 

/ [y, -y;i|x, / [y, -y/pds 

JO Jo 

I \Xs-K?ds. 

Jo 



We denote by Wg 



Ws := exp —{Hs, Ys) + $(s, y,')) 
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E 


snp\Yt-Y^\^ 


< E 


sup 




t<to 




t<to 



f Ws{Ys - y;, cj{s, Xs) - ct{s, Xi))dB, 
Jo 















2' 


+ E 


sup 


fws\Ys 






Ys))dBs 






t<to 


Jo 






















2' 


+ E 


sup 


fws\Ys 






Y',)dBs 






t<to 


Jo 











'■to 
+ / E 

} 

+ / E 



sup ly^ - Yi 

t<s 



/|4 



sup \Xt - X[ 

t<s 



/i4 



ds + 
ds 



to 



E 



sup\Yt-Y^\^\Xt-X't\^ 

t<s 



ds 



Then by noting that the three first terms are submartingales, we can apply Doob's inequality 
in order to obtain (with Cauchy- Schwartz inequality) for every to ^ ^■ 



E 



sup \Yt - Yl 

t<to 



/i4 



< E 



sup 

t<to 



[\ys - Y^,a{s,Xs) - a{s,X',))dB, 
Jo 



fto 

+ / E 



fto 

+ / E 



sup\Yt-Y^f 

.t<s 

snp\Xt-X^f 

t<s 



to 

ds+ I E 
ds 







snp\Yt-Y^\^\Xt-X't\^ 

t<s 



ds 



< E 



+ 



to 



Y,-Yf\X,-X'fds 



to 



E 



sup \Yt - Yl 

.t<s 



/i4 



+ 



to 



E 



sup \Xt - X[ 

t<s 



r\4 



ds + 
ds. 



to 



E 



sup\YtX\^\Xt-Xl\^ 

t<s 



ds 



The inequality xy < + implies 



E 



sup \Yt - Yl 
\t<to 



/|4 



to 



E 



sup|yt-y/| 

t<B 



ds + 



to 



E 



sup I Xt — X^ I 

t<s 



ds 



which with Gronwall's Lemma, gives us 





sup|yt-y/|4 


fto 






E 


< / E 


snp\Xt-X^\^ 


ds 




t<to 


Jo 


_t<s 





This concludes the proof of Proposition 137 
We finish this subsection by results without requiring the regularity of the set- valued map C 



□ 



Theorem 3.9. — Consider an admissible set-valued map C, varying in an absolutely continu- 
ous way. 

Let f,a:Ix Mf^ M*^ be bounded and Lipschitz with respect to the second variable: there exists 
a constant L such that for all t £ I and x,y £ 



and 



\a{t,x) - a{t,y)\ + \f{t,x) - f{t,y)\ < L\x - y\ 
\f{t,x)\ + \a{t,x)\ <L. 
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Then 121\) is well-posed. That means: for all initial data uq € C(0), there exists one and only 
one process {Xt)t£X solution of i21\) (in the sense of pathwise uniqueness). 

Proof: Theorem 13.91 is a consequence of Theorem 13.61 Indeed it suffices to repeat the proof 
developed in |31) . which permits to obtain Theorem 5.1 |31j (with abstract compactness ar- 
guments) as a consequence of Theorem 4.1 [31j (which corresponds to the well-posedness of 
deterministic Skorohod problem). □ 

Remark 3.10. — As already pointed out in Remark 5.1 |31| . the existence part of Theorem 
1 3. PI only requires the houndedness and the continuity of functions /, o". The Lipschitz regularity 
is necessary only for the uniqueness. 



3.2.2. Stability results. — Now we are looking for stability results: we let the stochastic pertur- 
bation a goes to and we prove that the corresponding solution tends to the solution of the 
deterministic sweeping process. 

Theorem 3.11. — Let ct'^ : Z x M'^ — > M'^ be maps satisfying there exists L with for all e > 

\cr^{t,x) - a%t,y)\ < L\x - y\ 

and 

II'^''IIl°°(IxR'*) — ^■ 

We assume that cr*^ tends to in L°°{I x W^) when e goes to 0. Under the assumptions of 
Theorem \3.8l we consider a fixed initial data uq G C{0). For each e G (0, 1], we denote (X|)tgj 
the (unique) process solution of 

( dXl + ^{C{t),Xl) 3 f{t,Xt)dt + a%t,Xl)dBt 

I Xl = uo. 

Denote x the solution of 

(IT 

-{t)+N{C{t),xit))3f{t,x{t)) 

x{0) = Uq , 

given by Theorem 1.1 of [18j and consider the deterministic process: defined for all (jJ ^ Vt by 
Xt{uj) = xt- 

Then X'' converges to X in L^(J7, L°°(X)).' 
for some constant Cuq independent on e. 

Proof: We denote F(e, •) the map defined on B := L'^in, L°°{I)) into B as follows: [F{e,X)]^ 
is the unique solution of the Skorohod Problem 

F{e,X)t+Kt = uo+ f a\s,X,)dBs+ ['f{s,Xs)ds. 

Jo Jo 

Analogously we define F(0, •). 

Then in the previous subsection, we have proved that X"" is the unique fixed point of F{e, •) and 
similarly X is the unique fixed point of -F(0, ■). In order to apply the "fixed point theorem with 
parameter", we check that the map F is continuous at e = 0. 

By Proposition [321 there exists a large enough integer p (only depending on the above constants 
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and L) such that for ah e > 0, the map F{e, is a ^-Lipschitzean map on B. We also deduce 
that 



< 



1 



< - X' -X „ + 
- 2 II lie ^ 



which gives 

\\X'-X\\^<2 F^\e,X)-F^P\0,X) 

It also remains to bound this quantity. By Proposition 13.21 the sequence {F{e,X))f: is uniformly 
bounded in B and the map F{e, •) are uniformly Lipschitz in B (with a constant denoted by k). 
Hence, Proposition 13.31 (see below) implies 

|F(2)(e,X) -F(2)(0,X)||^ 

< ||F(e, F(6, X)) - F{e, F(0, X))\\^ + \\F{e, F{0, X)) - F(0, F(0, X)) ||g 

< kk\\X\\B\\(T^\\Lco{XxR'i) + ^lk1lLcx,(2;xRd) 11-^(0, 

Then by iterating the reasoning, we deduce that there exists a constant Cuq (depending on uq 
through X, F{0,X), F^p\o,X) and on p) such that 1^ holds. □ 

Proposition 3.3. — With the notations of Theorem \3.11l there exists a constant k (indepen- 
dent on e) such that 

\\F{e,X) - F(0,X)||g < k\\XMa'\\Lo.(^i^^,y (32) 

Proof: Let us denote the processes := F{e,X), Y := F{0,X) and Z := F{e,X) - F{0,X). 
Then Z satisfies the following stochastic differential equation: 

dZt = - [dKl - dK^) + a'{t, Xt)dBt. 

With the same arguments as in the proof of Proposition 13. 2| it can be shown 



\Zt 



\Yt'-Yt\^< 



[\y,' -Ys,a'{s,Xs))dB, 
Jo 



+ 



— Y 1^(7* -I- II/t'^II^ 



Then applying Doob's inequality, we get for every tQ G I 

f{Y^ - Y,,a\s,Xs))dB, 
Jo 



E 


1 V |4 
sup \Zt\ 


< E 


sup 




t<to 




t<to 



+ 



to 



E 







sup iZtl"^ 



t<s 



/7 c -1- Wrr^ 11^ 

"T IF IIl°°{XxI 



< ll^'ll? 



+ 



L°°(IxI 
to 



^E 



to 



\ZJ ds 



E 







sup \Zt\ 

.t<s 



ds + llfj 



e||4 



lL°°{XxR'*)- 

By Cauchy-Schwartz inequality, we obtain with another constant c (depending on 



E 



sup \Zt\ 



t<to 



< 



_e||4 



lL°°(XxI 



+ 









/ E 


sup \Zt\ 


ds 


Jo 


_t<s 





By Gronwall's lemma, we obtain ([32 



□ 
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4. Euler scheme for stochastic sweeping process 

We refer the reader to the work of Y. Saisho [31], where the convergence of some discretized 
Skorohod problems to the continuous Skorohod problem is studied with a constant set C. For 
/ € C^{T) and a partition of X given by (t^ = ?T-/i)i<n<T/i-i) denote by Xh the following 
function: 

J Xh{t) := no, for ah t € [0,^] 

1 Xh{t) := Pc [xhK) + m - l{tl)] , for ah t € K,tl+\ 
Then in [31], it is proved that x^ strongly converges in L°°{I) to the unique solution x of the 
Skorohod problem 

f x(t) + kit) = lit), t G X 

{ (33) 
[ x(0) = uo. 

Moreover, the associated map := I — x^ strongly converges in L°°{I) to k. 

This section is devoted to the extension of such results for time-dependent sets C(-). Let us 
consider a set-valued map C(-) on Z = [0, T]. We now define a discretized solution Xh as follows: 

Xh{t) := no, for all t € [0,t\] 

Xh{t) := Pc(,n+i) [if^it-) + l{t) - l{tl)] , for ah t e]tl,tl+\ ^^^^ 

If the set-valued map C varies in an absolutely continuous way and takes uniformly prox-regular 
values, then it can be checked that Xh{t^) + l{t) — lif^) is close to C{t^'^^) and so its projection 
is single-valued for h small enough. The scheme is also well-defined for h small enough and 
semi-implicit as we consider C{tl+^). It is a prediction-correction algorithm: predicted point 
Xh{tl) + Kt) - Kth)^ that may not belong to C(tJJ+^), is projected onto C(t5J+^). Moreover, 
we refer the reader to the works (mentioned in the introduction) dealing with deterministic 
sweeping process. Such schemes are well-known in the framework of sweeping process (when 
I € W^'^il)) and corresponds to the so-called Catching-up Algorithm introduced by J.J Moreau 
in 



Remark ^.1. — If C{-) is a set-valued map taking uniformly prox-regular values and varying 
in an absolutely continuous way, it is well-known that the Euler scheme is convergent for smooth 
functions I (as in this case, Skorohod problem corresponds to sweeping process), see for example 
|18LI19| ." for I € W^'^{I), the sequence {xfi)h strongly converges to the unique solution x of (E^)- 

Theorem 4-2. — Consider an admissible set-valued map C, varying in an absolutely contin- 
uous way. For every function I G C^{I), the discretized solutions {xh)h strongly converges in 
L°°{I) to the unique solution x of 



Proof: According to Proposition 14. H {x)h is uniformly continuous with respect to h on Z. By 
Arzela-Ascoli's Theorem, {xh,h > 0} is relatively compact in C^{I). Then by uniqueness of 
solution to the continuous problem ([33|) . it suffices to show that every convergent subsequence 
(still denoted by {xhjh) converges to a solution of (p3]l . Let us fix such a convergent subsequence 
and write x for its limit. Obviously x is continuous. Since for t G]t^,t^"^^], x/i(t) € C{t^^), we 
deduce that 

dcit){xh{t)) < dH{C{tl+^),Cit)) < \v{tl+') - v{t)\ — > 0, 

where v is the continuous function given by Definition 13.41 So we get that for all t G I, 
x{t) € C{t). Moreover Proposition 14.11 implies that k = I — x has a bounded variation on X so 
it suffices to check (fT6|) . Indeed the algorithm implies that with kh ■= I — Xh and t (^]t^,t'^'^^] 

kh{t) - kh{tl) eN{C{tl+'),£h{t)). 
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We let the details to the reader and we refer to Theorem 4.1 of [31j for arguments in order 
to take the limit in the previous discrete inclusion and to deduce (fT6|l . The idea rests on the 
hypomonotonicity of the proximal normal cone. 

This permits to prove that x is solution of ([33]) and so by compactness we deduce that Xh 
converges to the unique solution of ([33|) . □ 

Remark 4- 3. — We would like to describe another way than the compactness argument, al- 
lowing us to prove the previous Theorem. Indeed by a similar reasoning than the one used in 
Proposition 14-11 we can obtain a discrete version of (6) Lemma 1.1 in [23j .• let 1,1 two contin- 
uous functions, we denote x and Xh, the continuous and discretized solutions of the Skorohod 
problem (SkP,l), then 

\\xh ~ ^hllioc(i) ^ 11^ ~ ^IIl°°(x) + IK ~ 

Moreover the implicit constant does not depend on h and depends only on I and I. So we conclude 
that the maps (/ x^) is Holder continuous of order ^ from (7*^(1, M"^) into itself on compact sets 
(and the bound can be chosen independent on h). Thanks to Proposition \3.1\ the map (/ — )• x) 
is Holder continuous of order ^ from C^{X,W^) into itself on compact sets too. 
So for I G C^{Z) and e > 0, there exists a smooth function I such that \\l — 1\\l'^(x) < £• Then 
we have 

\\xh - x\\l^(x) < P/i - XhWh^ii) + ll^h - x\\l^{x) + \\x - 
^From the Holder regularity of the maps (/ — )• Xh) and (/ — )■ x), we know that the first and third 
terms are bounded by e^^"^ (up to a numerical constant). Thanks to Remark \4.1\ the second term 
tends to asl(z W'^'^iT). Thus we deduce that Xh uniformly converges to x. 

Proposition 4-i- — The continuous functions x^ are uniformly continuous on X and x^ — I 
has a bounded variation, both them uniformly with respect to h. These estimations only depend 
on I via its L°°-norm and its uniform continuity modulus. 

This proposition can be seen as a "discrete version" of ()19p (or Lemma 1.1 (7) and Lemma 1.2 
of |23) ). so we will follow their proofs. 

Proof: We take again the notations, introduced in the proof of Proposition 13.11 (second step). 
We set kh := I — x^ and let m be a fixed integer. 

First step: Boundedness of a discrete variation \kh\{Tm+i) — \kh\{Tm). 

For every s,t G [Tm,Tm+i], thanks to the admissibility property, we get (with no,ni satisfying 
tl° £[s,s + h[ and t^' £]t - h, t]): 

\xh{t) - Xh{s)\ + \l{t) - l{s)\ > {ihis) - Xh{t) + l{t) - l{s),UiJ 

jti — 1 

> E i-^hK^') + ^f^i^h) + m^') - miuij 

n=nQ 

+ {-Xh{t) + Xh(.tT),UiJ + {-Xh{tl°)+Xh{s),u,J 
+ {l{t)-l{tl^),u.J + {l{tl^)-l{s),u,J 

rai — 1 

>sY, \ih{tt') - xu{ti) - i{ti+') + im 

n=no 

+ 6\xh{t)-xh{tl')-l{t) + l{tl')\ 
+ 5\xUO-Xh{s)-l{tl<^) + l{s)\. 

We have used that for all integer n and t G]t^,t^'^^] 

XH{t) - xh{tl) - lit) + m) e -^^{C{tl+'),XH{t)). 
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So we conclude 

\kh\{t) - \kh\{s) < i {\ih{t) - Xh{s)\ + \l{t) - l{s)\) , 
where we set the discrete variation 

ni —1 

Ikhlit) - \kh\is) := \kh{t) - khitl')\ + \kh{tl+') - k^m + \kh{s) - khitl 

n=no 

Consequently, 

\kh\iT^+i) - \kh\iT^) < K, (35) 

for some numerical constant, as for Proposition 13.11 
Second step: Uniform continuity of Xh- 

Let ?7 be a constant of prox-regularity of all the sets C(-). Furthermore for s,t € [Tm,Tm+i] we 
write 

\xh{t) - x;,(s)|2e-[l'=''IW-l'='^l(^)l/^ = \xh{tl') - x;,(t;;»)|2e-[l^''l(C)-l^'.l(C)]/'7 + Rest. 

We denote S for the first term and we only study it as the Rest can be similarly estimated. By 
a discrete differentiation, it comes 

ni — 1 

n=no 
ni — 1 

^ |-^(^n) _^^(^no)|2 J^g-[|fc;,|(t;:+^)-|fc;,|{0]/r; _g-[|fch|(i^)-|fch|(C)]/'? 

n=no 

Furthermore 

ni — 1 

^ = E {xh{C')-XHK),Xh{C') + Xf,{t^^)-^^^^^ (36) 

n=no 
ni — 1 

n=rao 

:= Si + 52. (38) 
Then in Si, we make the following replacement in the first term in the inner product 

x.(tr^) - x.iti) = {i{tt') - m)) - {k,{q+') - k,{ti)) . (39) 

Thus we write Si = I + II. The first term / equals to 

ni — 1 
n=no 

- mt^) - lK),h{tl+') + khitl) - 2A;,(t;^«))e-[l'=''l(tr^)-|fc.l(C)]/'y. 

The first quantity Ii (by producing the reverse manipulations with / instead of Xh) can be 
bounded by sups<j^<t2<t I ^(^2) — ^(^i)P by (|35]l . With the help of a change of variable, it can be 
shown that 

ni — 1 

i2<J2 - w)\ \m) - /(C)| < imt) - ms)] sup - iiti)\. 

n=no s<ti<t2<t 

We also have estimated the first term /, it remains to deal with the second one // (due to ()39p ). 
We recall that 
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The hypomonotonicity property of the proximal normal cone (see Definition 12.51 and Proposition 
EH) yields 

Indeed, the point G Q{t^°) and possibly does not belong to Q{t^^^) so we have to replace 

it by its projection onto which makes appear the two last quantities. We produce a 

similar reasoning for x/j(t^), in noting that 

We deduce that the second term is bounded by 

ni-l 
^ n=no 

+ c [\khm - |A:,|(.)] ([\vm - \v\{s)] + [\v\{t) - \v\{s)f + sup - l{h)\ 

\ S<ti<t2<t 

where c is a numerical constant and v is given by Definition 13.41 {v controls the variation of the 
set C7(.)). 

It remains to study the term S2 corresponding to ([37]) . Using 



^-[\kh\(t':)-\kh\(t?0)]/rj ^-[\kh\(t'!;^'^)-\kh\(C)]/v 

\kk\{tl^')-mq) 

it comes 
II + S2 

< [\khm - \kH\{s)] ([\vm - \v\{s)] + [\v\{t) - \v\{s)f + sup m) - . m 

\ S<ti<t2<t / 

Finally with (|35p . we obtain that 
\xh{t) - 

< sup \l{t2)-l{h)\^+ sup \l{t2)-l{h)\ + [\v\{t)-\v\{s)] + [\v\{t)-\v\{s)]\ (41) 

s<ti<t2<t s<ti<t2<t 

So we get 

sup \Xh{t) - Xh{s)\^ < sup \l{t2) - l{tl)\^ + sup \l{t2) - l{tl)\ 

Tm<S<t<Tm + l S<tl<t2<t S<tl<t2<t 

+ [\v\{Tm+l) - \v\{Tm)] + [\v\{Tm+ 

As in Proposition 13.11 the uniform continuity of v implies that [|?;|(Tm+i) — |f|(rm)] can be 
assumed small with respect to r (it suffices to take r sufficiently small). Hence, the collection 
of indices m is finite and so (1351) becomes 



|A;;,|(r)-|A;^|(0) <K' (42) 

for some numerical constant K' depending only on / (by its uniform mudulus continuity and 
its L°°-norm). Dividing the time-interval I with subintervals of length r, we deduce that the 
total variation \kh\{T) is (uniformly with respect to h) bounded by a constant depending on 
/ G C70(Z). 

Then ()42l implies that ()4ip holds for every s,t £ I, which gives the uniform continuity of Xh- D 
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Remark 4-4- — have detailed the bounds of the terms I, II and S2 in order to make 
explicit the extra term in ^0\ ). In the continuous versions of these results with time-independent 
set C{t) = C, it is well-known that we get II + S2 < (see (7) in Lemma 1.1 of |23j for 
example). 

This new term makes appear the "variation^' of the set C{-) in the corresponding interval [s,t\. 



Then we deduce the following result: 

Theorem 4-5. — Consider an admissible set-valued map C, varying in an absolutely continu- 
ous way and /, a be functions satisfying the assumptions of Theorem \3.9l Let X be the unique 
process, solution of 

( dXt + N(C7(t), Xt) 3 fit, Xt)dt + Xt)dBt 

\ Xq = Uq, 

where uq £ C(0) is a non- stochastic initial data, given by Theorem \ 3.9l For all w G we 
construct the following discretized process: 

- for all t G [0, tl] 

X^{oj):=uo, for all t e [0,tl] 

-for allte]tl,tl+'] 

X^{u:) :=Pc(i;j+i) [ifn M + /^/(t;:, X| M) + cT(t^, X| H) (S, - i?,.) 
Then for almost every w G il, we have: 



lim 



X^{oj)-X{uj) 



0. 



We refer the reader to Theorem 5.1 of [31j (mainly Lemmas 5.1 and 5.2 of this work) for a 
detailed proof of such result, using compactness arguments. 



5. Example of applications with particular set-valued maps 

In this section, we deal with some particular moving sets C, defined as an intersection of 
complements of convex sets. More precisely, we check that the assumptions (made in the previous 
theorems) are satisfied in the framework of |39j . Let us recall it. 

We consider the Euclidean space M'^, equipped with its euclidean metric | |, its inner product 
(•, •) and B the closed unit ball in M'^. Let I := [0,T] be a bounded closed time-interval and 
for i G {1, ..,p} let gi : X x M'^ — )■ M be functions (which can be thought as "constraints"). We 
introduce the sets Qi{t) for every t £ I by: 

Qi{t) := {xGM^ g^{t,x) > o} , 

and the following one 

p 

1=1 

which represents the set of "feasible configurations x" . We suppose that for all t G X and all 
i G {1, ..,p}, gi{t,-) is a convex function. We suppose there exists c > and for all t in [0,T] 
open sets Ui{t) D Qi{t) verifying 

dH(Qi(t),M'^\[/i(t)) >c, (AO) 
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where dn denotes the Hausdorff distance. 

Moreover we assume that there exist a,/3,M,K > such that gi G (X x {Qi + kE)) and 
satisfies: 

Vi G X, X G Ui{t), a < \V^gi{t,x)\ < (3, (Al) 

ytel, xeUiit), \dtgi{t,x)\ < ^, {A2) 

and 

Vt G X, X G Ui{t), \Dlgi{t,x)\ < M. {A3) 

Vt G X, X G Ui{t), \dtV^gi{t,x)\ < M. {A4) 

We denote by 

I{t,x) := {i, gi{t,x) = 0} 

the set of "active contraints" and 

Ip{t,x) := {i, gi{t,x) < p} , 

for some p > 0. As explained before, we want to deal with admissible sets Q{t) and so we have 
to make the following important assumption: there exist constants p, 7 > such that for all 
t G X and all x G Q{t) 



(Rp) 



Xi\S/xgiit,x)\ < 7 

i£lp{t,x) 

for every nonnegative coefficients A,. We refer to |39| for a first use of this kind of "reverse 
triangle inequality" (i?o) and {Rp). 

We want to apply the previous results to the set-valued map Q{-) in order to get well-posedness 
results for the following stochastic sweeping process 

r dXt + N{Q{t),Xt)3 f{t,Xt)dt + a{t,Xt)dBt 

{ (43) 
i Xo = uo£ Q(0). 

By Proposition 2.8 in [39], we know that for all t G X, Q{t) is uniformly prox-regular and we 
can describe its proximal normal cone. 

Proposition 5.1. — Under the assumption (Rq), there exists a constant rj such that for all 
t £ I, the set Q{t) is rj-prox-regular. 
Moreover, for allt^X and x G Q{t), 

N(g(i),x)= Yl N(Qi(i),x) = - Yl ^'^^x9i{t,x). 

i€l{t,x) i€l{t,x) 

In order to study the set-valued map Q(-), we need this technical lemma 

Lemma 5.1. — There exist constants v,T,p\r > such that for allt and x G Q{t): there 
exists u satisfying: 

- \u\ = l 

— for aZ/ s G [t — r, t + r] n X, y G -B(x, 2r) and i G y), 

{Vxgi{s,y),u) >iy. (44) 

Remark 5.2. — Such result was already proved in Lemma 2.10 [39j with r = r = 0. 
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Proof: Let t € X and x € Q{t), we set the fohowing cone 

Np(Q(t),x) :=- Yl ^^'^x9^{t,x) 

i£lp(t,x) 

and its polar cone 

Cp{Qit),x) := Np(Q(t), x)° := {w; G G Np(Q(t), x) , (t;, u;) < o} . 

According to the classical orthogonal decomposition of a Hilbert space as the sum of mutually 
polar cones (see [28j ). we have: 

= PNp{Q(t),x) + PCp{Q{t),x), 

where P denotes the Euclidean projection. 

So let us consider for i G Ip{t,x) the corresponding decomposition of Vx giit,x): 

Vx gi{t, x) = ai + bi£ Np(Q(t), x) + Cp{Q{t),x). 

Assumption (|Aip gives us: \bi\ < \Vx gi{t,x)\ < 13. Since G l>Sp{Q{t),x), it can be written: 
= — ^ AjVq gj{t, x) , Aj > involving 



l^il = \'^x9i{t,x) -ai\= ^ XjVxQjit.x) + (1 + \i)Vxgi{t,x) 
Then using the inverse triangle inequality and Assumption (jAip . we get: 



a 



As a consequence, it comes: 



a 



< \bi\ < /3. 



(45) 



Since G Np{Q{t),x) and = PNp{Q(t),x)(Va; 5i(*, a;)), we obtain 

2(6i, -Va;5i(t,x)) = \bi - Vxgi{t,x)\^ - \bi\^ - \V ^ gi{t,x)\'^ 
= |oiP - l^iP - |V^5i(t,x)p 



I, i9 

< - 6,- P < ^. 



7- 



Now we set 



u := 



ieip{t,x) 



eCp{Q{t),x). 



This is well-defined because ()46p and Assumption (jylip imply that for any j G Ip{t,x) 



i&Ip(t,x) 



/3 



\i€lp{t,x) I 

Then (jM]) for (s,?/) = (t,x) follows from (gS]) and (gH]) with 

v2 



a 
2^ 



(46) 
(47) 

(48) 



a 



For /)' = p/2, Lipschitz regularity (jAip and (|^2p imply for y G B{x, 2r) and s G — r, t + r] flZ 
with r,r < p/(8/3) 

Ip'{s,y) C Ip{t,x). 

Thanks to Assumptions ^A2>\f and (041), it can be shown 

(V^ gi{s, y),v)>u' - M{2r + r) > i/, 
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with V = v' 12 and r,r < v/{2M). 

Consequently (glD holds for p' = p/2, u = and r, r < min{p/(8/3), z^/(2M)}. □ 

Proposition 5.2. — For all t = [0,T], i/ie set-valued map Q{-) is admissible. 

Proof: The uniform prox-regularity of Q{t) is asserted in Proposition 15.11 It also remains to 
check Property p^ . 

Let t (z I and x G Q{t) and r, r, n given by the previous Lemma. We fix s E [t — T,t + t] OZ 
and y S -B(x, 2r) n dQ{s). Every u G N((5(s), y) fl B can be written as follows 

v = - ^ K^xQiis.y) 

i£l{s,y) 

with nonnegative coefficients Aj. By Assumption (j^ip . 

/3 5] A, >L 

iG-r(s,y) 

Then from Assumptions l\A3\i and (031) ; we deduce a Lipschitz regularity for the gradient V^; gi 
and so, we get 

iel{s,y) 
iel{s,y) 

> v/P. (49) 

We set 5 := u/ f5 and Ux ■= —u. 

So we have proved that for every x G dQ(t) we can find Ux, satisfying for all y G dQ{s)r\B{x, 2r) 
with |t - s[ < T and u G N(Q(s), y) n B 

(w,tx^)>'5. (50) 
Then Property (|15p is obtained by choosing a bounded covering of \J dQ{s) C M'^ with 

S, \s~t\<T 

balls of radius r. □ 

Proposition 5.3. — The set-valued map Q is Lipschitz continuous with respect to the Haus- 
dorff distance. 

This result was already proved in [39]. We give a proof for an easy reference. 
Proof: Consider t,s G [0,T] and x G Q{t), let us construct a point close to x belonging to 
Q{s). Let u given by Lemma l5.H we introduce z{h) := x + hu with /i > 0. We claim that for 
h < hi := min(c, p/{j3 + v)) , 

Vi G {1, gi{t,z{h)) >hv 

(where c and v are introduced in ()AOp and (j44p ). Indeed for h < c, due to the convexity of 
gi(t,-), it comes 

gi{t,x + hu) > gi{t,x) + h{Vxgiit,x), u). 
As a consequence, for i G Ip{t,x), 

gi{t, X + /in) > hi', 
by ([33]) ■ Moreover for every i ^ Ip{t,xq), according to Assumption l\Al\i . 

gi{t, xq + hu) > p — hp > hv 
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Figure 1. The sets (3(0) and Q{i) for some i > 0. 



\ih < . Thus for h < hi, we have gi(t,x + hu) > hi' for all i € {1, --jp}- That is why we 

p + 

deduce from Assumption iA2^ that 2:(/i) € Q{s) if /ii^ > /3|t — s|. Setting £ := , if |i — s| < i, 

it can be written that 



(3 

dQ(s){x) < inf \x — z{h)\ < —\t — s\. 

hh'>l3\t—s\ V 



Consequently we obtain if |t — s| < 



/3, 



c^H(Q(i),<5(s)) = max sup sup dQ(t)(x) <-|t-s|. 

\a;GQ(t) x(^Q(s) ) ^ 

This inequahty is actually satisfied for any t, s € [0,T]. To check it, it suffices to divide the 
corresponding interval into subintervals of length I and to apply the triangle inequality. □ 

Remark 5.3. — We would like to show the importance of Assumption (Rq) for the Lipschitz 
regularity of Q with an example. Consider the two-dimensional space M? (d = 2) and chose two 
constraints 

gi{t,x):=X2 and g2{t,x) := —e~^^ — X2 — t. 

The two constraints are smooth and convex and satisfy the assumptions Ml|) -i H4)) . Then it is 
easy to check that for allt >0 the line {x2 = 0} belongs to the set Q{t) (see FigureU^. However 
for t > 0, we have 

Q{t) C {{xi,X2), X2 >0 and xi < — log(t)} . 
So we deduce that dH{Q{t),Q{0)) = oo for all t > 0. Computing the gradients, we get 

V.^gi{t,x) = (0,1) and V,^ g2{t, x) = {-e~''\ -I). 

Assumption (Rq) is not satisfied as we have: 

\V,,gi{t,x) + V^g2{t,x)\ = e'^^ 

which could be as small as we want. 



Conclusion: In this framework (Assumptions ()AOP - ()^4p and {Rp)), we have proved that the 
set-valued map Q is admissible and Lipschitz continuous. 

We can also apply Theorem 13.91 and obtain well-posedness results concerning the following 
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stochastic differential inclusion (with maps / and a which are bounded and Lipschitz with 
respect to the second variable): 

' dXt + N(g(t), Xi) B fit, Xt)dt + Xt)dBt 

< Xt(^Q{t) 

^ Xo = Uq, 

for all Uq G Q(0). Moreover, an Euler scheme can be used to approach the solution, as it is 
convergent thanks to Theorem 14.51 

This numerical scheme ()34p makes appear the Euclidean projection Pq(j) onto the set Q{t). As 
the set Q{t) is not assumed to be convex but only uniformly prox-regular, we have not efficient 
numerical algorithm in order to compute this projection. 

In the case of the set Q, defined (see the beginning of this section) as the intersection of com- 
plements of convex sets, an implementable scheme has been proposed in |24|, 125^ I26| and then 
extended to moving sets Q{t) in |39j . The idea is to replace Q{t) with a convex set Q(t,x) 
(depending on the variable x) defined for any point x G Q{t) as follows: 

Q{t,x) :={yGM^ gi{t,x) + {Vg^{t,x),y - x) >0 Vi} . 

This set is convex and included into Q{t), due to the convexity of the functions gi{t,-). This 
substitution is convenient because classical methods can be employed to compute the projection 
onto a convex set. Yet this replacement raises some difficulties for the numerical analysis which 
are solved in the framework of deterministic sweeping process, in proving that Q{t,x) is a good 
local approximation of Q{t) around the point x (we refer the reader to j39j for a detailed 
study and corresponding results). The new algorithm (instead of (|34p) is the following one: for 
I G C^{I) we now define a discretized solution Xh as follows: 

f ihit) ■■= no, for all t G [0,tj^] 

\ Ut) ■■= PQitl-^\,,itl)) i^hiq) + m - m)] , for all t G]q,tl+']. ^^^^ 

It would be interesting to prove its convergence. In |39j . the second author has already proved 
its convergence in the framework of sweeping process for / G W^'^(Z). The extension of this 
result to continuous functions / G C^iT) is still open. Indeed, the main difficulty rests on the 
lack of regularity of the map x — )• Q{t,x) (even for time-independent constraints gi). 

We finish this work by briefly presenting a stochastic model of crowd motion, which is an 
extension of a deterministic one introduced in [27j . We refer the reader to [38\ \2fi\ 12 7| for 
a complete and detailed description of this model, which takes into account the direct conflict 
between people. 

We consider N persons identified to rigid disks. The center of the i-th disk is denoted by G 
and its radius by r^. Since overlapping is forbidden, the vector of positions q = (qi, .., qA?) G 
has to belong to the "set of feasible configurations" , defined by 

Q := {q G M^A^, A,(q) > Vi / j} , (52) 

where -Djj(q) = |qi — Qjl — (?"?. + rj) is the signed distance between disks i and j. 

By denoting by U(q) = {Ui{qi), ..,Un{c[n)) G I^^^ the global spontaneous velocity of the 
individuals, the crowd motion model can be written 

dq + N(Q,q) 3 \J{q)dt. 
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Now people's hesitation or panic can be modelled by a stochastic perturbation a{t, ci)dBt, which 
gives the following stochastic model: 

dq + N(Q, q) 3 U(q)dt + ait, ci)dBt, (53) 

associated to a Brownian motion (i?i)t>o- 

In order to get well-posedness results, it also suffices to check that Q is admissible. This was 
already proved in [32j by a direct approach and later in [39] via showing Assumption (Rp) (since 
Assumptions ()AOp . (j^dip and ^A3^ can be easily satisfied). 

Theorem 5.4- — If\J,(T{t,-) are Lipschitzean and a is bounded, then i53\) is well-posed: for 
every initial condition q(0) € Q, there is a unique solution to i53\) . 

Remark 5.5. — Since we deal with a constant set Q, [31\ already allows to conclude. However 
the current results allow to get the well-posedness of i53\) when the radii depend on time under 
the following assumptions: Vi are uniformly time-Lipschitz and 

inf inf n(t) > 0. 

te[0,T] i 
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